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c/3 The excesses of the cosmic positron fraction recently measured by PAMELA and 
the electron spectra by ATIC, PPB-BETS, Fermi and H.E.S.S. indicate the existence 

j> ■ of primary electron and positron sources. The possible explanations include dark 

<n ; 

(NJ , matter annihilation, decay, and astrophysical origin, like pulsars. In this work we 

in 

show that these three scenarios can all explain the experimental results of the cosmic 
(N ■ + 

e excess. However, it may be difficult to discriminate these different scenarios by the 

oo : 

. local measurements of electrons and positrons. We propose possible discriminations 

among these scenarios through the synchrotron and inverse Compton radiation of 

X' 

' the primary electrons/positrons from the region close to the Galactic center. Taking 



typical configurations, we find the three scenarios predict quite different spectra 
and skymaps of the synchrotron and inverse Compton radiation, though there are 
relatively large uncertainties. The most prominent differences come from the energy 
band 10 4 ~ 10 9 MHz for synchrotron emission and > 10 GeV for inverse Compton 
emission. It might be able to discriminate at least the annihilating dark matter 
scenario from the other two given the high precision synchrotron and diffuse 7-ray 
skymaps in the future. 

PACS numbers: 95.35.+d,95.85.Ry,95.85.Pw,97.60.Gb 
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I. INTRODUCTION 



The anti-matter particles such as positrons and antiprotons in cosmic rays (CRs) are 
very important in understanding the origin and propagation of CRs. These particles are 
usually produced by the CR nuclei that interact with the interstellar medium (ISM) when 
propagating in the interstellar space. Precise measurements of these particles will provide 
us valuable information about the primary CR sources and the interaction with matter. In 
addition, these secondary anti-matter particles have relative lower fluxes and characteristic 
spectra. Therefore they make themselves good objects to study the exotic origin of CRs, 
such as from dark matter (DM) annihilation or decay. 

Recently the PAMELA collaboration released their first CR measurements on the positron 
fraction 1| and p/p ratio 2[. The positron fraction of PAMELA data shows an uprise above 
~ 10 GeV till to ~ 100 GeV, which exceeds the background estimation of the conventional 
CR propagation model. This result is consistent with previous measurements by, e.g., HEAT 
[s| and AMS Q]. On the other hand, the p/p ratio is compatible with the background. The 
results of PAMELA strongly indicate the existence of such sources that mainly generate 
leptonic particles at this energy range. The electron spectrum up to several TeV reported by 
ATIC collaboration also shows an obvious excess with an interesting bump around 300 ~ 800 
GeV |5fl . In addition, the measurements of the electron spectrum bvPPB-BETS M, H.E.S.S. 

in Q 

3, 18|], and most recently by Fermi |9[ all show the excesses of electron spectra, although they 
are not fully consistent between each other. The studies show that the PAMELA result on 
the positron fraction is consistent with the electron spectrum measured by ATIC or Fermi 
assuming the equal amount of the production of positrons and electrons (e.g., [id. 111!]). 

One possible primary electron/positron source in the Galaxy is the pulsar. Pulsars and 
their nebulae are well known cosmic particle accelerators. Although the quantitative details 
of the acceleration processes are still open for study, early radio observations have established 



them as important high energy cosmic electron and positron sources (e.g., 12j). X-ray and 
7-ray observations show that some of the accelerated particles can reach an energy of a few 
tens of TeV, and there are indications that the particle distribution cuts off in the TeV energy 
3.g., [yl, 14, 15|). These particles produce emission over a broadband frequency range 



range (e. 

in the source region, which has been observed and studied extensively (e.g., 16[). Some of 
these particles will escape into the ISM becoming high energy cosmic electrons and positrons. 
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It is shown that one or several nearby pulsars will be able to contribute enough positrons to 



reproduce the PAMELA data 



17 



18) as well as ATIC data (l9| 



Another potential candidate, which has been widely discussed to solve the 



positron/electron excess problem, is the DM annihilation (e.g., 20J]) or decay (e.g., 21 



23|) 



It has been shown in many studies that the local measurements of the electrons/positrons 
can not provide enough power to discriminate the DM origin and the pulsar origin of the 
electrons/positrons. On one hand, the ATIC data suggest a sharp cutoff near ~ 600 GeV. 
Such a feature can naturally appear for a DM origin when closing to the mass of DM particle, 
if the DM annihilates or decays directly to e + e~ pair, as shown in Ref. However, the 
same thing is also expected if the observed excess is dominated by a single pulsar and/or its 
nebula and the accelerated high energy electrons and positrons have experienced significant 
energy loss as they propagate in the source and in the ISM [lj|. The narrow peak resulting 



from local pulsars is even indistinguishable from that of DM origin [24]. On the other 
hand, if the electron spectrum is as smooth as that from Fermi, there are also degeneracies 
between the pulsar and DM interpretations. A continuously distributed pulsar population 
in the Galaxy might be able to generate smoother electron spectrum due to different energy 
losses of various pulsars. While for the DM scenario, if the electrons and positrons come 
from DM annihilation or decay, the spectrum can also be smooth, and is consistent with the 



Fermi measurements 11 ] 



In this work we study the photon emission associated with the different scenarios to see 
how the degeneracies between these models can be broken. Our point is that although the 
different scenarios can give almost the same local electrons/positrons, the spatial extrapola- 
tions of them might be significantly different, especially in the Galactic center (GC) region. 
The pulsars mainly concentrate in the Galactic plane, while the DM is spherically distributed 
in an extended halo and highly concentrated in the center of the halo. In addition, DM an- 
nihilation is proportional to p 2 (r), different from that of DM decay (oc p(r)). Therefore we 
anticipate the existence of remarkable differences in the skymaps of the synchrotron and 
inverse Compton (IC) emission from these extra primary electrons/positrons for the three 
scenarios. 

We compare the differences of synchrotron as well as IC radiation among the pulsars, 
annihilating DM and decaying DM models, on the premise that they all reproduce the 
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positron fraction and total electron + positron flux data. The spectral and directional 
features of the radiation are discussed in detail. We compare the results with the WMAP 
haze data 25, 26] and EGRET diffuse 7-ray data j^], and give the observable signals by 



future experiments such as Fermi/GLAST. 

This paper is organized as follows. In Sec. II we will briefly present typical model con- 
figurations of the three scenarios which can fit the data simultaneously. The discussions 
on the synchrotron and IC emission are given in Sec. Ill and IV respectively. In Sec. V 
we discuss the uncertainties of our predictions. Finally, we give a summary and draw the 
conclusion in Sec. VI. 

II. DIFFERENT SCENARIOS TO ACCOUNT FOR THE e ± DATA 



We start with the propagation equation of electrons/positrons in the Galaxy 

uih o o 

= Q(r,p) + V ■ - Vc^) + ^-/ D vvq--^ 

-^[#-|(V-V c ^) 



(1) 



where ip is the number density of CR particles per unit momentum interval, Q(r,p) is 
the source term, describing the primary particles injected into the interstellar medium, 
D xx is the spatial diffusion coefficient, V c is the convection velocity. The reacceleration 
process is described by the diffusion in momentum space, with the diffusion coefficient 



D 



- 3( 5(4_ ( 52)(4_ ( 5) wZ) — , where va is the Alfven speed, w is the ratio of magnetohydrodynamic 
wave energy density to the magnetic field energy density, which characterizes the level of 
turbulence, p = dp/dt is the momentum loss rate, mainly induced by synchrotron radiation 



and IC scattering for the electrons/positrons for energies > 10 GeV 



In this work, 



we use the numerical package GALPROP 



28j to calculate the propagation of the primary 



electrons/positrons and the background CRs. 

We have adopted the conventional propagation model where all the CR data are re- 
produced by the model. The cosmic nuclei secondary to primary ratio, such as B/C, the 
unstable secondary to stable secondary, such as 10 Be/ 9 Be, and the local proton and electron 
spectra are taken to constrain the propagation parameters. A conventional diffusion + con- 
vection (DC) propagation model [22| is adopted. For details of the propagation model see the 



discussion in Ref. 



22| . The propagation parameters are taken as: half height of the propa- 
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gation halo z^ — 4 kpc; diffusion coefficient D xx = j3D (p/ p ) s with D = 2.5 x 10 28 cm 2 s" 1 , 
5 = 0.55 for rigidity p > 4 GV and 5 = for p < 4 GV; the convection velocity is adopted 
as a linear function of coordinate z with V^(z = 0) = and dV c /dz = 6 km s _1 kpc -1 . The 
reacceleration is not included. 

The source term Q(r,p) is different for various scenarios. For DM annihilation, the source 
function of electrons/positrons has the form 

(av} A p 2 (r) dN(E) 



Qa{ y i E) = BF- 



2 m DM 



dE 



(2) 



where BF represents the "boost factor" which can come from the dumpiness of DM sub- 



structures (e.g., 



29 



30 



3l| ) and/or the so-called Sommerfeld effect 20[. In the former 



case, BF is energy and spatial dependent, BF ~ BF(r, E). However, the study based on 
N-body simulations shows that the DM substructures do not tend to strongly enhance the 
DM annihilation signals 32] . Thus in this work we adopt the Sommerfeld effect as the 



"boost factor" for discussion. (<jv)a is the velocity weighted annihilation cross section be- 



dN(E) 



A IS 



fore being enhanced by the Sommerfeld effect, ttibm is the mass of DM particle, 
the electron/positron spectrum from one pair of DM annihilation, and pir) is DM density 
distribution in the Galaxy. To satisfy the relic density constraint we adopt the cross section 
(ctv)a = 3 x 10 26 cm 3 s . For the DM density profile we adopt the Merritt profile 33] 



p(r) = p exp 



2 

a 



Da 



(3) 



with a = 0.2, r„ 2 = 25 kpc and the local DM density p & = 0.3 GeV cm -3 [34| at r = R Q = 
8.5 kpc. 

For DM decay, the source term is given by 

1 p(r) dN 



Qd{t,E) 



TDMm DM dE 



(4) 



D 



where tdm is the life time of DM, ^ \d is the electron/positron spectrum from the decay 
of one DM particle. 

The source term for Galactic pulsars is 



dN 

Q P (R,z,E) = K-f(R,z)- — 



(5) 



where K is the normalization factor representing the total luminosity of Galactic pulsars, 



f(R,z) is the pulsar spatial distribution and -j]?\p is the average electron/positron energy 
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spectrum radiated from the pulsars. The spatial distribution of pulsars can be parameterized 

as 



f(R,z)<x(^- 



exp 



b(R - R & 
R Q 



(6) 



where R is the Galactocentric radius, and z is the distance away from the Galactic Plane. 
Different from spherically symmetric form of DM distribution, pulsars are mainly concen- 
trated at the Galactic Plane with z s ~ 0.2 kpc. The primary electron/positron spectrum 
injected by pulsars is generally assumed to be a power law form with an exponential cutoff 
at high energies 

dN 

— =E~ a exp(-E/E cut ), (7) 




FIG. 1: Left: The positron fractions predicted by the three scenarios after solar modulation, 



compared with HEAT [3|, 



ELA [jj data. Ri 



H.E.S.S. 



and 



the three scenarios, compared with observations of ATIC 5[, PPB-BETS 
Fermi 9j. We assume DM annihilates or decays into pure lepton final states with equal branching 
ratios to e + e~, /i + fi~ and t + t~. 



We have taken appropriate parameters in each scenario to give best fit to the PAMELA 
and ATIC or Fermi data. For the annihilating DM, we adopt ttidm = 1 TeV, and a boost 
factor BF = 800. The annihilation channels are assumed to be pure leptons with equal 
branching ratios to e + e~, fi + fi~ and t + t~ respectively. We find that the PAMELA and 
ATIC data can not be fitted well by only one lepton flavor final state. Note, however, for 
the case of Fermi data the pure or r + r~ channel might be able to fit the data. We will 

further discuss this point in Sec. V. For the decaying DM, we set DM mass ttidm = 2 TeV 
and life time r ~ 1.08 x 10 26 s. The decay channels are the same as the annihilation scenario. 
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For pulsars, we adopt the spatial distribution parameters a = 1.0, b = 1.8 35j, and the 
power law index of the energy spectrum a ~ 1.4 with an exponential cutoff above ~ 800 
GeV. The spectrum index 1.4 can give good description to the Fermi data. However, to 
better fit the ATIC data we need to introduce harder spectrum and/or super-exponential 
cutoff. 

The positron fraction and the total electron + positron flux are shown in Figure [H The 



solar modulation effect is calculated using the force field approximation [37J with a potential 
of 500 MV in this work. For the background electron choice, we adopt the DC model given 
in Ref.[22]. A rescale factor 0.9 on the electron normalization is adopted to better fit the 
PAMELA and ATIC/Fermi data simultaneously. This model is actually similar with the so- 



called conventional model [38|, which fits all the pre- Fermi CR data and is also successfully 
used to model the diffuse gamma-ray spectrum measured by Fermi at intermediate Galactic 
latitudes [39]. 



Finally we should point out that we take continuous distribution of the sources (pulsar or 
DM) to recover the PAMELA and ATIC/Fermi data in this work. Note that it is also possible 

nn 

that the locally observed electrons/positrons come from one or several nearby pulsars |17U19l] 



or DM clumps 



401 ] . For the pulsar scenario, this assumption does not affect the discussion 



since the extrapolation to the whole Galaxy will hold regardless of how much the far pulsars 
can contribute to the local electron flux. For the annihilating DM model, although a single 
DM clump which is close enough to the Earth may be possible to give large boost factor, it 
is found to be of little probability to survive in a realistic DM distribution model 4jJ, |42 |. 
For the decaying DM scenario, since the "boost" effect from DM decay will be much weaker 
we will expect an even smaller probability to find a clump close and massive enough to 
explain the electron/electron data, compared with the annihilating DM scenario. Therefore 
the assumption of continuous distribution of sources is reasonable. 



III. SYNCHROTRON RADIATION FROM THE THREE SCENARIOS 

In this section we study the synchrotron emission from the primary electrons and positrons 
in the three scenarios to account for the locally observed cosmic e excesses. The predictions 
are based on the same configurations as discussed in the last section. 

The synchrotron emissivity with frequency v of a single electron with Lorentz factor 7 in 



s 



a magnetic field B, is given as 43j 
e(>,7) = 4:V3nr e m e cu L x 2 ^K i/3 (x)K 1/3 (x) - y [#4/3 fa) ~ K i/z( x )] J er S s s~ 1 Hz" 1 (8) 

where r e = e 2 /m e c 2 is the classical electron radius, x = vji^rfvi) with vl = eB /2irm e c 
the Lamor frequency, and Ki(x) is the modified Bessel function of order /. For an ensemble 



of electrons/positrons with power law spectrum oc E~ a , the corresponding synchrotron 
intensity is also a power law form l[y 
strength in the Galaxy is adopted as 



oc v 13 , with index f3 = (a — l)/2. The magnetic field 



44] 



R — Rq\ ( \z\ 



B(R,z) = B exp f exp ( -L± \ , (9) 

where Rb = 10 kpc and zb = 2 kpc, and the local magnetic field strength B = 5 /xG. Such 
a magnetic field can match the 408 MHz synchrotron longitude and latitude distributions 



Solving Eq. ([T]) we get the electron/positron distributions in the Galaxy. Then we can 
derive the synchrotron emissivity at any point in the Galaxy The line of sight integral of the 
emissivity gives the synchrotron flux as a function of direction (/, b), which can be compared 
with the observation. 



A. Results from the three scenarios 



synchrotron spectra: GC _ synchrotron spectra: AntiGC 




10 3 10 4 10 5 1 6 1 7 1 B 10 9 10 10 10 4 10 5 10 6 1 7 1 a 10 9 



v, MHz v, MHz 



FIG. 2: The average synchrotron spectra of the three scenarios within a bin size of 20° x 20° around 
GC and anti-GC, respectively. 

In Figure [21 we give the average synchrotron spectra of the three different scenarios in a 
20° x 20° window centered at the Galactic center (GC) and anti-GC respectively. In the GC 
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direction the annihilating DM scenario has the largest signal, while in the anti-GC direction 
the decaying DM model has the largest signal though not very high in absolute fluxes. It 
is shown that these models have quite different synchrotron radiation spectra in the GC 
direction from 10 4 MHz to 10 9 MHz, even though they have similar contributions to the 
local electrons/positrons. 

We give the results by smoothing the signals within the 20° x 20° window centered at 
the GC so that we can diminish the uncertainties of our prediction. The largest uncertainty 
comes from the DM density profiles. The prediction of DM annihilation at the GC from 
cuspy or cored density profile may be orders of difference. By smoothing the signal within 
a large window the uncertainty is considerably decreased (see detailed discussion in Sec. 5). 
Further with a large window we can treat the pulsars as continuous distribution with average 
emission spectrum. Uncertainties from background are also deceased in a large window. 




a 




latitude, degree 



180160140 120100 80 60 40 20 340 320 300 280 260 240 220 200180 

longitude, degree 



FIG. 3: The synchrotron latitude profile for |/| < 10° (left) and longitude profile for |6| < 10° 
(right), at the frequency of 61 GHz. The solid black line indicates the contribution from background 
electrons and positrons, and the other lines represent the total synchrotron intensities from the 
three scenarios together with the background respectively. 



Given the large difference of synchrotron radiation from 10 4 MHz to 10 9 MHz in Figure 
[21 we choose a frequency randomly in this range, 61 GHz, to demonstrate the synchrotron 
longitude and latitude profiles in Figure [3j It is shown that the DM model will have larger 
gradient around GC, especially for the longitude profile, when compared with the pulsar 
model. For annihilating DM model the longitude profile is extremely cuspy around / = 0. 
This is a significant feature to distinguish among these models. In addition, the absolute 
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fluxes of the synchrotron radiation in the inner Galaxy are also different. For pulsar model, 
the longitude and latitude profiles almost follow the background distributions and might 
be difficult to be distinguished from the background. It is not hard to understand these 
profile features. The pulsar distribution is similar to the source distribution of CRs such 
as the supernova remnants, while the DM distribution is spherically symmetric and much 
more concentrated at GC than the pulsars. In addition the DM annihilation products are 
proportional to the square of DM density, p 2 (r), while the decay products depend only on 
the DM density p(r). If the DM density is cuspy in the center of the halo, like that shown 
by many numerical simulations, the annihilation scenario will show the most cuspy profile 
at the GC direction. 



B. The WMAP haze 

The WMAP satellite has made all-sky survey in the microwave bands with several fre- 
quencies in the range ~ 10—100 GHz. The Galactic foreground dominates the radio 
emissions of these skymaps. Such emissions are expected to be produced by standard in- 
terstellar medium emission mechanisms, such as thermal dust, spinning dust, ionized gas 
and synchrotron radiation. After excluding all known contributions, however, Finkbeiner 
et al. reported an anomalous excess of the microwave emission around the GC which is 



called WMAP "haze" 25| . In Ref . 26] , the authors confirmed the existence of the haze by 



analyzing WMAP three-year data and gave the possible implication of its origin. It is also 



45]. 



proposed that the haze might come from the DM annihilation 

In Figure H] we compare our calculated synchrotron profiles in the inner Galaxy from 
the three scenarios with the WMAP data. From this figure we can see that annihilating 
DM model may reproduce the WMAP haze data by taking proper density profile, while the 
decaying DM and pulsar models can not produce the steep rise toward the GC as shown in 
data. 



IV. DIFFUSE 7-RAYS FROM IC SCATTERING 



In this section we study the diffuse 7-ray emission in the same three scenarios as discussed 
in the last two sections. Diffuse 7-rays can be produced through IC scattering of energetic 
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FIG. 4: Residual synchrotron intensities after subtracting the background for the three scenarios, 
compared with WMAP Haze data [26I ] at five different frequencies. 



electrons and positrons with the interstellar radiation field (ISRF). The ISRF consists of 
radio emission of cosmic microwave background (CMB), infrared emission from dust and 
optical emission from stars. A model to calculate the ISRF distribution based on the realistic 



Galactic stellar and dust distributions was built by Porter and Strong 



46j | and this result 



has been included in the GALPROP package. In this work we use GALPROP to calculate 
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the IC- induced diffuse 7-rays generated by the primary electrons/positrons as well as the 



background diffuse 7-rays generated by conventional CRs |44(. Note that the final state 
gamma ray radiation of annihilating or decaying DM by internal bremsstrahlung, which is 
significant near the mass of DM, is not taken into account here. As shown in Ref. [l(| the 
IC component will dominate the diffuse 7-rays in the Fermi energy range, which is also of 
major interest in this work. 

In Figure O we show the diffuse 7-ray spectra for two sky regions, |/| < 30°, \b\ < 5° and 



< 60°, \b\ < 10° for the three scenarios, corresponding to the region A 38| and region H 



27l | when discussing the EGRET data. The IC 7-rays from DM annihilation are significantly 



larger than the other two scenarios, especially for energies higher than several GeV. The 
contributions from pulsars are very small and almost undetectable if the uncertainties of the 
background are taken into account. The signal of decaying DM model is moderate. The 
results from all the three scenarios do not violate the observations of EGRET. 

For comparison we also plot in Figure [5] the contributions of annihilating DM with 
Navarro- Frenk- White (NFW) profile, 

P(r) = /' 12 , (10) 



1 + 



r 



with r s = 20 kpc and the local DM density p & = 0.3 GeV cm -3 at r = R Q = 8.5 kpc. 
Since NFW profile is steeper in the halo center than Merritt profile, it will give harder 
7-ray spectra. In a smaller sky region, the differences between these two kinds of DM 
distribution profiles are more obvious, as shown in Figure [6j We can see that the IC 7-rays 
and synchrotron radiation can provide a good test of the annihilating DM scenario and even 
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give constraints on the inner profile of DM distribution 

Figure [7] gives the diffuse 7-ray longitude and latitude profiles in the inner Galaxy, for 
two energy intervals, 4 — 10 GeV and 10 — 100 GeV respectively. The data points in the top 
two panels are from EGRET diffuse skymaps. It is shown that the precise measurements of 
the 7-ray gradient near GC will be also very helpful to discriminate different models, and 
possibly reveal the particle properties of DM {^J. 
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FIG. 5: The contributions to the diffuse 7-ray spectra in the inner Galaxy of the three scenarios 
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FIG. 6: Same as Figure [5] but for annihilation DM within a smaller region |/| < 5°, |6| < 5°. NFW 
and Merritt profiles are adopted to show the difference. 

V. DISCUSSION ON THE UNCERTAINTIES OF PREDICTION 



In the last two sections, we study the prospects to discriminate the three scenarios ac- 
counting for the PAMELA and ATIC data, i.e. annihilating DM, decaying DM and Galactic 
pulsars, through the synchrotron and IC radiation. For each scenario, a specific model is 
chosen. In this section we will discuss the uncertainties of the model configurations. 
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FIG. 7: The diffuse 7— ray latitude (average with |/| < 30°) and longitude (average with |6| < 5°) 
profiles. Top two panels are for energies 4 — 10 GeV, and the bottom two panels for 10 — 100 GeV. 



A. Particle physics model of DM 

In this work we adopt the annihilation/decay modes of DM only to leptons with equal 
branching ratios to e + e~, fi + fi~ and r + r~ to explain the PAMELA and ATIC data. However, 
our predictions of the synchrotron and IC radiation are not generally based on this specific 
choice. Since no matter how a model is constructed, the final electron and positron spectra 
have to be adjusted to fit the PAMELA and ATIC data. It is the electron and positron 
spectra that determine our prediction of the synchrotron and IC radiation. 

Actually for a flavor blind DM the 1:1:1 branching ratios to e + e~, and t + t~ 

is also the most natural and typical choice. First the PAMELA antiproton data severely 
constrained the gauge boson and quark final states 2j. As the first order approximation we 
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FIG. 8: The positron fractions and electron spectra in the case of pure e + e , or t + t final 

state of DM. 

assume pure leptonic final states. Further, in Figure Owe show that the PAMELA and ATIC 
data are not fitted very well if only e + e~, or t + t~ final state is assumed. If we gauge 

the positron ratio to the PAMELA data, the electron/positron spectrum will be too hard to 
account for the ATIC data if only taking the e + e~ channel, while it is too soft for the t + t~ 
channel. In addition, pure or t + t~ could not reproduce the sharp falling at ~ 600 

GeV at the ATIC spectrum. A combination of the three channels is necessary to account for 
the data from the two experiments simultaneously. The e + e~ channel is mainly to produce 
the sharp falling at ATIC data while the r + r~ channel contributes to the PAMELA data at 
low energy. 

According to the latest published data by Fermi Q and HESS 8], the bump feature of 
electron spectrum is not as prominent as that observed by ATIC. In this case, the direct 
channel to e + e~ to describe the fast drop of electron spectrum is not necessary. A DM model 



with purely as the final states has been proposed to well fit the new data [ll|. While 

from Fig. [S]we can see that a pure r + r~ channel might also be able to fit the Fermi/H.E.S.S. 
data. In the upper panels of Fig. [9] we show the positron fractions and total (e + +e~) spectra 
from DM annihilation and decay, of which the final states have equal branching ratios into 



fi + fi~ and t + t~ 1 , such as in the gauged L^ — L T model 50J. In this calculation we adopt the 
same NFW DM density profile as in EqflU} m^M = 1.5 TeV for annihilation and tbdm = 3 
TeV for decay. The predicted synchrotron and diffuse 7-ray spectra from the GC region are 
shown in the bottom panels of Figure We can see from this figure that the annihilation 



1 A pure or t + t channel can yield similar results, with the photon emission almost unchanged. 
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FIG. 9: The predictions of the positron fraction (upper-left), total (e + +e~) spectrum (upper-right), 
synchrotron spectrum (bottom-left) and diffuse 7-ray spectrum (bottom-right) for a 1.5 TeV DM 
annihilation model and a 3 TeV DM decay model, both of which with equal branching ratios into 
and t + t~. The pulsar scenario is the same as in Figures [2] and 

DM model could still be discriminated. 

It should be noted that the current measurements of ATIC and Fermi are not fully 
consistent. Fermi result has high precision, however, both are dominated by systematic 
errors. The Fermi data seem to favor an astrophysical origin, such as from nearby pulsars. 
However, the ATIC data favor the DM scenario more. From our study we notice that the 
DM annihilation scenario seems always be distinguishable from the other two. 



B. Spatial distribution of DM 

The key point of the present work is that the three scenarios have quite different spatial 
distributions, especially near the Galactic center, although all the three scenarios are eas- 
ily to account for the local measurements by PAMELA and ATIC. It may be difficult to 
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discriminate the three scenarios by local observations near the Solar system. However, the 
difference of the three scenarios near the GC should be large if we can try to observe their 
products, such as the synchrotron or IC radiation as we studied in this work. Therefore, the 
most important uncertainties should come from the spatial distribution of DM density. 

In order to remove the uncertainties of the DM density profile we have averaged the 
synchrotron and IC radiation over a large window around the Galactic center. This will 
smooth the signals around the GC and diminish the uncertainties from DM density profile 
efficiently as well as the distribution of pulsars. 

Numerical simulations generally suggest cuspy DM profiles in the center of the DM halo 



[e.g. 



51 



53]). In this work we adopt the Merritt profile with intermediate central 



density between the cuspy profiles like NFW and the cored one like isothermal. This DM 
profile gives a good description of the WMAP haze data, as shown in Figure HI 

To show the uncertainty of DM profile here we consider the most conservative case with 
the cored isothermal profile. The density profile is given as 

with a = 5 kpc and the local DM density p & = 0.3 GeV cm -3 at r = R & = 8.5 kpc. With 
the same local density p Q , the cored isothermal profile is the least cuspied distribution. Thus 
we anticipate it will produce the weakest photon signals at the GC. We show in Figure [TU] 
the synchrotron spectra within a bin size of 20° x 20° around the GC and the diffuse 7-ray 
spectra around the GC for cored isothermal DM profile compared with the pulsar model and 
the background. It is shown the differences between these models indeed become smaller 
compared with that given in Figure [2j The results shown in Figure [10] could be taken as 
the lower limit of our prediction, if other related models and parameters are fixed. Any 
other DM density profile will lead to a more clear discrimination of these predictions. In 
this conservative case with cored isothermal profile, it might be not easy to discriminate the 
DM signals. 



C. Propagation model 

Another source of uncertainties of our prediction comes from the propagation models. 
Actually requiring the propagation model to reproduce all the CR data gives strong con- 
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FIG. 10: The same as Figure 3 but for the isothermal DM profile. 

straints on the propagation parameters. The secondary to primary ratio, such as B/C, the 
unstable secondary to stable secondary, such as 10 Be/ 9 Be, and the local proton, electron 
data are used to determine these parameters (see the discussion of the propagation model in 



Ref 221]) . Under these requirement our results given in the previous sections will not change 
much by adjusting the propagation parameters, since the propagation models give similar 
descriptions to all kinds of CRs. Once other species satisfy data the model should also give 
similar electron and positron distributions at the same time. We will explicitly show the 
expectation in the following. 

In the previous sections we adopt the conventional model that all the cosmic ray spectra 
are consistent with observations. We have ignored the effect of reacceleration in Eq. and 
tept the effect of convection, which leads to better fit to the PAMELA data at low energies 



22j | . We will show in the following the propagation parameters will modify our results very 
slightly. 

As an example we give the synchrotron spectra for diffusion + reacceleration (DR) prop- 



agation model adopting the same propagation parameters as in [22j, which, contrary to 
the propagation model adopted in the present work, keeps the effect of reacceleration while 
ignores the effect of convection. In Figure [TT] we plot the synchrotron spectra in the DR 
propagation model for the three scenarios to account for the PAMELA and ATIC data. We 
see that these results are almost the same as the prediction in the DC model shown in the 
left panel of Figure [2j 

Another important uncertainty of propagation model comes from the height of the dif- 
fusion region. We have studied a propagation model with diffusion height 2kpc. The prop- 
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FIG. 11: The same as Figure 3 for Left panel: DR propagation model and Right panel: DC 
propagation model but with 2kpc diffuse halo height. 

agation parameters are also adjusted to reproduce all the cosmic ray observations. In the 
right panel of Fig. [11] we show the synchrotron predictions in this 2kpc diffuse halo height 
propagation model. We see that the change of propagation parameters generally affect our 
predictions slightly. 



D. Other astrophysical uncertainties 

The synchrotron and IC radiation directly depend on the magnetic field and ISRF respec- 
tively. The GALPROP model has adopted the latest result of the magnetic field and ISRF, 
which are given according to the most recent astronomical observations. To fully figure 
out the uncertainties of these astronomical observations is not easy. However, taking these 
uncertainties into account will not change our conclusion in this work, since the change will 
affect all the three scenarios simultaneously. For example the strength of the magnetic field 
may have the uncertainty of a factor 2. To reduce the magnetic field twice will reduce all the 
synchrotron radiation of the three scenarios, together with the background, twice. This is 
equivalent to say that all the curves in the synchrotron plots should shift downwards, while 
the relative differences between these models are not changed. Similarly change of ISRF will 
shift the IC radiation for all the three scenarios at the same time. However, observations of 
radio or diffuse 7-rays by EGRET and Fermi especially at high latitude will constrain the 
magnetic field or ISRF that can not be too small. 
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VI. SUMMARY AND CONCLUSION 

To explain the recent measurements of the positron fraction by PAMELA and electron 
spectrum by ATIC, Fermi and HESS, several scenarios including the DM and pulsars are 
proposed. In this paper we studied the perspective to discriminate these models using the 
synchrotron and IC radiation generated by these electrons/positrons. The point is that 
although various kinds of models degenerate in the local environment 2 and give similar 
electrons/positrons spectra, their different spatial extension might be revealed by the ac- 
companied photon emission. 

By properly adjusting the model parameters we first build three benchmark models which 
can reproduce the measured electron/positron data. We then studied the synchrotron and 
IC radiation from these primary electrons and positrons. We find that around the GC region, 
the synchrotron spectra for frequencies from 10 4 ~ 10 9 MHz and the IC 7-ray spectra from 
several GeV to several hundred GeV are significantly different among the three scenarios, 
especially between annihilating DM scenario and others. For the directional profiles, the 
DM models show larger gradient near the GC, especially in the longitude profiles, than the 
pulsar model and the background. The photon emissions for pulsar model are very small 
and similar with the background, and will be the most difficult one to be detected. While 
for DM models, especially the annihilating DM model, there are strong observable signals. 
Finally we discuss the possible uncertainties of these conclusions. The major uncertainty 
comes from the DM inner profile. For the most conservative case with cored isothermal DM 
profile, the signals of DM scenarios become weaker and might be hard to be discriminated. 
However, the isothermal profile is generally disfavored by numerical simulations. Other 
uncertainties, such as the DM annihilation or decay channels, the propagation parameters 
and other astrophysical inputs, seems not to change the relative differences among the signals 
of the three scenarios significantly. 

It should be pointed out that in this work the boost factor for the annihilating DM 



scenario is assumed to be universal, such as the Sommerfeld effect. 



the boost effect is spatially dependent, e.g., from the DM subhalos [29|, |30j, l31|- If so, the 



t is also possible that 



enhancement near the GC region will not be as important as at large radius in the halo. 



2 High energy electrons lose energy fast and the detected electrons should indeed be "local", e.g., within 
~ 1 kpc [54f | . 
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Actually, the calculation based on the N-body simulation indicates that the boost factor 
from DM substructure for anti-matter particles from DM annihilation in the Galaxy is 
generally negligible 32J, |55[ . A recent calculation on the DM dumpiness boost factor taking 
into account the Sommerfeld effect still finds no remarkable enhancement effect, except for 



fine tuning to the strongly resonant case 
generally held. 



Therefore the conclusions in this work are 
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